Effective radii of deuteron induced reactions 



Shintaro Hashimoto,' - * Masanobu Yahiro,^ Kazuyuki Ogata,^ Kosho Minomo,^ and Satoshi Chiba' 

' Advanced Science Research Center, Japan Atomic Energy Agency, 
Ibaraki 319-1195, Japan 
^Department of Physics, Kyushu University, Fukuoka 812-8581, Japan 
(Dated: January 19, 2013) 

The continuum-discretized coupled-channels method (CDCC) for exclusive reactions and the eikonal re- 
action theory (ERT) as an extension of CDCC to inclusive reactions are applied to deuteron induced reac- 
tions. The CDCC result reproduces experimental data on the reaction cross section for d+''*Ni scattering at 
200 MeV/nucleon and ERT does data on the neutron-stripping cross section for inclusive ^Li(ci, n) reaction 
at 40 MeV. For deuteron induced reactions at 200 MeV/nucleon, target-dependence of the reaction, elastic- 
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plained by simple formulae. Accuracy of the Glauber model is also investigated. 
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I. INTRODUCTION 

Understanding of the fusion reaction mechanism is one 
of the most important and challenging subjects in nuclear 
physics. The fusion reaction consists of complete and incom- 
plete fusion processes. In the complete fusion process, all of 
the projectile is absorbed by the target nucleus. In the incom- 
plete fusion process, meanwhile, a part of the projectile is ab- 
sorbed, while other part(s) of the projectile is emitted. The 
complete fusion process at low incident energies is essential 
to understand the production of superheavy nuclei. The in- 
complete fusion process in the scattering of unstable nuclei at 
intermediate energies is important to extract information on 
the projectile from the scattering. Actually, the nucleon re- 
moval reaction widely used for the spectroscopy of unstable 
nuclei [1] is composed of the nucleon-stripping reaction as a 
consequence of the incomplete fusion process and the elastic- 
breakup reaction as a result of the direct reaction process. Fur- 
thermore, the proton stripping process in inclusive ^Li((i, n) 
reaction at incident energies up to 50 MeV attracts wide inter- 
ests of not only nuclear physicists but also nuclear engineers, 
because emitted neutrons through the process are planned to 
be used in the international fusion materials irradiation facil- 
ity (IFMIF) [2]. Accurate evaluation of the proton-stripping 
cross section is highly required. 

The theoretical tool of analyzing the incomplete fusion 
process at intermediate energies is the Glauber model [3]. 
The theoretical foundation of the model is investigated in 
Ref. [4]. The Glauber model is based on the eikonal and 
the adiabatic approximation; the latter is known to make the 
elastic -breakup and removal cross sections diverge when the 
Coulomb interaction is included; see for example Ref. [5]. 
The Glauber model has thus been applied only for lighter tar- 
gets in which the Coulomb interaction is negligible; see for 
example Refs. [1, 6-10]. Very recently, inclusive ^Li((i, ?i) 
reaction at 40 MeV [11] was analyzed by the hybrid calcu- 
lation [12] in which the elastic -breakup component is evalu- 
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ated by the continuum-discretized coupled-channels method 
(CDCC) [13, 14] and the proton stripping component is by 
the Glauber model. The analysis was successful in reproduc- 
ing the data [11], even if the Coulomb interaction is neglected 
in the Glauber-model calculation. Such a hybrid calculation 
should be justified by more accurate reaction theories. 

CDCC is an accurate method for treating exclusive reac- 
tions such as elastic scattering and elastic-breakup reactions. 
The theoretical foundation of CDCC is shown in Refs. [15- 
17]. CDCC succeeded in reproducing data on the scattering 
of stable and unstable projectiles [13, 14, 18-31]. Very re- 
cently, CDCC was extended to inclusive reactions such as nu- 
cleon stripping reactions [32]. This method is referred to as 
the eikonal reaction theory (ERT). In ERT, the adiabatic ap- 
proximation is not made for the Coulomb interaction, so that 
the elastic -breakup and the nucleon removal reaction never 
diverge. ERT is thus applicable for both light and heavy tar- 
gets. ERT assumes the eikonal approximation to be good. The 
formulation starts with the eikonal approximation, but non- 
eikonal corrections are made by calculating fusion cross sec- 
tions with CDCC. This is essential progress in the theory on 
fusion reactions. 

Extensive measurements of total-reaction and nucleon- 
removal cross sections are now being made for the scat- 
tering of unstable nuclei at intermediate energies, say 100- 
300 MeV/nucleon, in MSU, RIKEN, and GSI. Accurate un- 
derstanding of the fusion reaction mechanism is thus highly 
required at intermediate energies. In this paper, we mainly 
consider deuteron induced reactions at 200 MeV/nucleon as a 
typical case and analyze integrated cross sections of the reac- 
tions with CDCC and ERT. Deuteron is fragile just as unsta- 
ble nuclei and furthermore it has no ambiguity on the struc- 
ture. In this sense, deuteron is the most suitable projectile 
to understand the fusion reaction mechanism. We will show 
that CDCC reproduces experimental data on the reaction cross 
section for rf+^^Ni scattering at 200 MeV/nucleon and ERT 
does data on the neutron-stripping cross section for inclusive 
^Li((i, n) reaction at 40 MeV. Target-mass-number (A) de- 
pendence of the reaction, elastic-breakup, nucleon-removal, 
nucleon-stripping, incomplete- and complete-fusion cross 
sections for deuteron induced reactions at 200 MeV/nucleon 
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is clearly explained with simple formulae. Accuracy of the with the operator 
Glauber model will be investigated. 

ERT is recapitulated in Sec. II. Numerical results of CDCC 
and ERT are presented in Sec. III. Section IV is devoted to 
summary. 



O = 



JKZ 



(5) 



II. EIKONAL REACTION THEORY 

A. Three-body model 

Deuteron (d) is the system in which proton (p) and neutron 
(n) are weakly bound. It is thus natural to assume that scat- 
tering of d from target T is well described by the p + n+T 
three-body model. Actually, the model is successful in repro- 
ducing the experimental data on elastic scattering and breakup 
reactions of d [13, 14]. The model Hamiltonian is 



2p. 



Vl + U{rp,rr,) + h 



with 



(1) 



(2) 



where h = Tr + V{r) denotes the intrinsic Hamiltonian of d 
that consists of the kinetic-energy operator T,, and the interac- 
tion V. Furthermore, /i is the reduced mass between d and T, 

Up^'' (U^'') is the nuclear part of the proton (neutron) optical 

(c) 

potential and Up is the Coulomb interaction between p and 
T. The three-dimensional vector R = (6, Z) stands for the co- 
ordinate between d and T, while r is the coordinate between p 
and n. The vector = {b^, z^-) for x = p or nis the coordi- 
nate between x and T. The total wave function '^{R, r) of the 
three-body system is then obtained by solving the three-body 
Schrodinger equation 



[H -E]'^{R,r) 



0. 



(3) 



In the three-body model, transitions of the incident flux to 
the inelastic (target-excitation) channels are expressed by the 



imaginary parts of Uj 

r(N) 



(N) 



and t/f). 



The imaginary part of 



Up denotes the absorption of p by T, while the imaginary 

part of [/jI^' corresponds to the absorption of n by T. There- 
fore, the three-body model implicitly assumes that p and n are 
absorbed independently. 



B. Separation of S'-matrix 

We consider d scattering at intermediate energies, say 
200 MeV/nucleon. Since the eikonal approximation is con- 
sidered to be good for the scattering, the S'-matrix elements 
and several types of cross sections are described by ERT. Ac- 
curacy of the eikonal approximation is investigated later In 
this subsection we recapitulate ERT for deuteron scattering. 
In the eikonal approximation, the three-body wave function ^I* 
is assumed to be 



(4) 



where K = yj2^,{E — h)/h and v = hK / fi are wave- 
number and velocity operators of the motion of deuteron rel- 
ative to T, respectively. When Eq. (4) is inserted into Eq. (3), 
we have a term including V^i/) but it is neglected, since ip is 
slowly varying with R compared with O. The neglect leads 
Eq. (3) to 



(6) 



Regarding Z as "time" and solving Eq. (6) iteratively, we ob- 
tain the formal solution 



ip = exp 



-iV 



dz'd^'uo 



(7) 



where V is the "time" ordering operator Taking Z to cx) in 
Eq. (7), we get the S'-matrix operator 



exp 



-iV 



dZO'^UO 



(8) 



In the Glauber model, h is replaced by the ground-state en- 
ergy Co of d. This adiabatic approximation reduces O^UO 
and V in Eq. (8) to U/ (hvo) and 1, respectively, where vq is 
the velocity of d in the ground state relative to T. This is noth- 
ing but the 5-matrix in the Glauber model. Thus, ERT is an 
extension of the Glauber model. 

The operator UO describes the change in the motion of 
p and n in d during the collision. The change is small for the 
short-range nuclear interactions, Up^^ and U^\ while large 
for the long-range Coulomb interaction Up . Therefore, the 
adiabatic approximation that neglects this change is good for 
the nuclear interactions but not for the Coulomb interaction. 

More quantitative discussion can be made by considering 
the matrix element 



dZ((/);^|Ot[/o|0g) 



between the ground state 
eo and its continuum state 



^i{Ko'K)Ru 



_ J 

(9) 

/)o of d with the intrinsic energy 
with the intrinsic momentum 
hk and energy e, where HKq {hK) is the momentum of d 
in the ground (continuum) state relative to T, and Ru is the 
range of the interaction considered. As an example, let us 
consider the d+^'^^Pb scattering at 200 MeV/nucleon. The 
spectrum of the elastic breakup reaction, dcEB / de, has a peak 
around e = 10 MeV. The interaction range Rjj is about 



while 



7.1 fm for the nuclear interactions, uf^'' and u'ii^\ 

(c) 

that for the Coulomb interaction Up is infinity. Hence, 
A = (A'o — K)Ru ~ 0.55 < tt for the nuclear interactions, 
but oo for Up . Since the adiabatic approximation is good 
for A ^ TT, the approximation is acceptable for the nuclear 
interactions, but not for the Coulomb interaction. Actually, 
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the breakup cross section is known to diverge in the adiabatic 
approximation [5, 22]. 

The fact that the adiabatic approximation is fairly good for 
U^'^ indicates that C/n^' is commutable with O. Therefore, 
we can make the replacement 

Ot{/(N)OoC/fV(^^^o). (10) 

The accuracy of Eq. (10) is confirmed later by numerical cal- 
culations. Using Eq. (10), we get 



S — SnSp 



with 



Sn = exp 



Sp = cxp 



-iV 
-iV 



(11) 



(12) 



(13) 



Thus, iS* can be separated into the neutron and proton parts, 5„ 
and Sp, respectively. The neutron part 5„ describes scattering 
of n by U^'' and recoil of p by the scattering. However, a 
velocity caused by the recoil is much smaller than the initial 
velocity vq of p, so that the recoil effect is negligible. Similar 
interpretation is possible for Sp. The operator Sp is the formal 
solution of the Schrodinger equation 



2^ 



h + Uj^Hrp) + UP{rp)-E 



*p = 0. (14) 



and Sn is the solution of the Schrodinger equation 



*n = 0. 



(15) 



Hence, the matrix elements of Sn and Sp can be obtained 
by solving Eqs. (14) and (15) with eikonal-CDCC [22] in 
which the eikonal approximation is made in CDCC calcula- 
tions. Non-eikonal corrections to 5„ and Sp can be easily 
made by using CDCC instead of eikonal-CDCC in solving 
Eqs. (14) and (15). 



C. Integrated cross sections 

In this subsection, we derive several formulae of integrated 
cross sections with the product form (11), following the for- 
mulation of the cross sections in the Glauber model [7, 33]. 
The reaction and elastic -breakup cross sections, ctr and cteb, 
respectively, are defined by 



d^b[l - \{^0\SpSn\^0)? 



(16) 



aEB = J d^b[{cbo\\SpSn\^\M-\{MSpSn\<l>o)\%(lV 

The cross sections ctr and cteb can be evaluated from the 
asymptotic form of ^I* that is obtained by solving Eq. (3) with 
CDCC. 



The total fusion cross section cttf is defined by 

tXTF = fTR - CTEB = J d^b{^o\[{l-\SpSn\^mo). (18) 

The total fusion cross section can be decomposed into the 
neutron-stripping cross section cr„:STR, the proton-stripping 
cross section (Tp:STR and the complete-fusion cross section 
CTcp: 



fTF = CTn-.STR + CpiSTR + ^CF, 



(19) 



where 

fT„:STR = I d^b{(j>o\\Spf{l-\SnfMo), (20) 

fTp:STR = f d'b{cf,o\\Sn\\l~\Sp\^Mo), (21) 
j2i,/ / i/i ic |2\/i Ic |2\ 



acF = J d'b{(t>o\il - |5„n(l ~ ISpl'Mo). (22) 

The factor 15^1^(1 — l^np) in cr„:gTR, shows thatp is scattered 
by T while n is absorbed by T, and the factor (1 — — 
jS'pp) in (JcF means that both p and n are absorbed by T. The 
sum of (JniSTR and (Tp:STR describes the incomplete fusion 
cross section ajp: 



ClF — CTniSTR + CTp-.STR 



(23) 



In the neutron removal reaction, n is either absorbed or 
scattered by T, while p is independently scattered by T. Hence, 
the cross section a^n is the sum of cteb and (Tn-.sTR- 

= Cteb + CT,i:STR- (24) 
The neutron-stripping cross section ct„:str is rewritten into 

CT„:STR - / d^biMiil - \SpSn\^) - (1 - ISpl^mo) 



where 



Cttf — cttf(p), 



o-tf(p) = ctr(p)-cteb(p) 



(25) 



(26) 



with 



Mp) = J d^b[l ~ \{(j>o\Sp\<j>o)\% (27) 
cteb(p) = f d^b[{<j>o\\Sp\^\M ~ \{MSp\M\% (28) 



Here, cttf(p), ctr(p) and cteb(p) are the total fusion, reaction 
and elastic-breakup cross sections induced by Up^^ + Up'^'^ 
only. The cross sections, ctr (p) and ctebIp), can be evaluated 
from the asymptotic form of 'i'p that are obtained by solv- 
ing Eq. (14) with CDCC. Thus, we can evaluate ct„ str with 
Eq. (25). 

Similarly, the proton removal cross section CT_p is obtained 

by 



CT-p — Cteb + CT„:STR 



(29) 
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and the proton-stripping cross section (Tp:STR is rewritten into 

CpiSTR = cttf - cttfI"-), (30) 

where 

(Ttf(?^) = crYi{n) - aEB{n) (31) 

with 

an{n) ^ J d'b[l - \{^o\Sn\M\% (32) 

^ebH = I d'b[{cbo\\Sn\''\(l>o) -\{cbo\Sn\4>o)\^]. (33) 

The cross sections, crR(n) and uebI"-), can be evaluated 
from the asymptotic form of that are obtained by solv- 
ing Eq. (15) with CDCC. We can thus evaluate (Tp:STR with 
Eq. (30). Finally, a^p is obtained from cttf and aip — 
CTniSTR + cTpiSTR by using the relation ctcf = <^tf - ctif- 

D. Tests of the eikonal and the adiabatic approximation 

In ERT, non-eikonal corrections to the integrated cross 
sections are taken into account by using CDCC in- 
stead of eikonal-CDCC. For the d+^Be scattering at 
200 MeV/nucleon, the correction is found to be less than 1 % 
for ctr, (Teb, cttiiStr, and (Jp.sTR- For the d+^^^Pb scattering 
at 200 MeV/nucleon, the correction is 1.5% for crpj, cteb and 
CniSTR and 16% for CTp:STR. Thus, the non-eikonal correc- 
tion is small except for (7p:STR for heavy targets. As shown in 
Eq. (30), CTp:STR is approximately obtained by 

cTpiSTR ~ ctr - o-R(n-), (34) 

since cteb ^ ctr - The 1.5% correction appears in fTR because 
of the strong Coulomb field, while the correction is negligible 
in CTR(n) as a consequence of the absence of the Coulomb 
field. Thus, one can conclude that the 16% correction required 
for crp:STR is nothing but the 1.5% correction for ctr. Note 
that (Jp.sTR is much smaller than ctr- Meanwhile, ct„:str is 
given by 

a'„:STR ~ ctr - cTRb) (35) 

The 1.5 % corrections appear in both ctr and <Jn{p). The can- 
cellation between the two corrections makes the non-eikonal 
correction small for (T„:str- 

In ERT, the adiabatic approximation is assumed to be good 
for the nuclear potential uj^\ This can be tested by setting 

U{R, r„) = UPiR) + [/f + U^fKr„) 

in the Schrodinger equation (3). In this setup, the projectile 
breakup is induced only by U^\rn), since the argument Vp 

of C/p^^ and t/p^'' has been replaced by R. Switching the 
adiabatic approximation on the Schrodinger equation corre- 
sponds to the replacement (10). For the d+^Be scattering 
at 200 MeV/nucleon, the error due to the approximation is 
0.3% for (TR and 2% for cteb- For the d+^°^Ph scattering at 
200 MeV/nucleon, the error is 0.4% for ctr and 6% for <teb- 
Errors due to these approximations are even smaller for heav- 
ier projectiles such as ■^^Ne [32]. 



III. NUMERICAL RESULTS 

We use the Koning-Delaroche global optical potential [34] 
as [/p^'' and and the Ohmura potential [35] as V that 

reproduces the deuteron binding energy |£o| = 2.23 MeV 
As the model space of CDCC calculation, s-, p-, and d- 
wave breakup states with k < 1.0 fm~^ are taken. Each k- 
continuum is divided into small bins with a common width 
Afc = 0.1 fm~^, and the breakup states within each bin are 
averaged into a single state. Maximum values of r and R are 
''max = 200 fm and i?,nax = 200 fm, respectively. 

A. d+^^Ni elastic scattering at 400 MeV 

In this subsection, we consider the d+^^Ni elastic scatter- 
ing at 200 MeV/nucleon, because the elastic cross section was 
measured and the reaction cross section was evaluated with 
the optical potential analysis [36]. 

Figure 1 shows the elastic cross section as a function of the 
center-of-mass (cm.) scattering angle 9. The solid (dashed) 
line represents a result of CDCC calculation with (without) 
the spin-orbit interactions of J/p^"* and t/,!^' . The solid line 
well reproduces the experimental data. Large deviation of the 
dashed line from the solid line for 6 > 10° shows that the 
spin-orbit interactions yield a significant effect on the elastic 
cross section. 




10 20 30 40 

eideg.] 

FIG. 1: (Color online) The elastic cross section of d+^'^Ni scattering 
at 200 MeV/nucleon as a function of the cm. scattering angle 6. 
The solid (dashed) line stands for result of CDCC calculation with 
(without) the spin-orbit interactions of the proton and neutron optical 
potentials. The experimental data are taken from Ref. [36]. 



In the Glauber model, the spin-orbit interactions and the 
Coulomb breakup are neglected, and furthermore the eikonal 
and the adiabatic approximation are made. The Coulomb 
breakup can be neglected by replacing Up'^\rp) by Up'^^R) 
in CDCC calculation. In Fig. 2, the dotted line is the result 
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of CDCC calculation with neglecting both the spin-orbit in- 
teractions and the Coulomb breakup. The dotted line agrees 
with the dashed line of Fig. 1, that is, the result of CDCC 
calculation with the Coulomb breakup and without spin-orbit 
interactions. Thus, the Coulomb breakup effect is small. In 
Fig. 2, the dot-dashed line represents a result of the Glauber- 
model calculation. The large deviation of the dot-dashed line 
from the dotted line comes from the eikonal and adiabatic 
approximations, more precisely from the eikonal approxima- 
tion. Eventually, the result of the Glauber model (the dot- 
dashed line) largely deviates from the full-CDCC result (the 
solid line) in which both the spin-orbit interactions and the 
Coulomb breakup are taken into account. Thus, the Glauber 
model does not work well for the elastic cross section for 
> 10°. 



breakup, since the Coulomb breakup becomes more signifi- 
cant for heavier targets. 
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TABLE I: Effects of the spin-orbit interaction and the Coulomb 
breakup on ctr, (teb and cttf for the d+^**Ni scattering at 
200 MeV/nucleon. The experimental data is taken from Ref. [36]. 
The cross sections are shown in units of mb. 
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FIG. 2: (Color online) Comparison of the Glauber model with 
CDCC for the elastic cross section of rf+^*Ni scattering at 
200 MeV/nucleon. The solid line is the same as the solid line in 
Fig. 1. The dotted line stands for the result of CDCC calculation 
with neglecting both spin-orbit interactions and Coulomb breakup. 
The dot-dashed line represents the result of the Glauber-model cal- 
culation. The experimental data are taken from Ref. [36]. 



B. Inclusive ^Li(d, n) reaction at 40 MeV 

The double differential cross section (DDX) of inclusive 
^Li((i, n) reaction was measured at 40 MeV [11]. The main 
part of the DDX consists of the elastic-breakup and the pro- 
ton stripping parts. When the elastic-breakup DDX is cal- 
culated with CDCC and subtracted from the measured DDX, 
the angular and the energy dependence of the remaining DDX 
is well reproduced by the Glauber model [12]. This indi- 
cates that the proton-stripping cross section (TpiSTR can be 
obtained by fitting the theoretical DDX calculated by the 
Glauber model to the remaining DDX, and integrating it 
over the angle and the energy. The cTpiSTR thus extracted is 
244 ± 34(tlicor.) ± 37(cxp.) mb; the theoretical error comes 
from ambiguity of the fitting. ERT gives o'p:STR = 253 mb 
and the Glauber model does 214 mb. Thus, the two theoretical 
results are consistent with the experimental data. The ERT re- 
sult seems to be slightly better than the Glauber-model result 
for this case. 



Relation between reaction and elastic-breakup cross 
sections 



The reaction cross section calculated by full-CDCC is 
1056 mb, while the value extracted from the measured elastic 
cross section with an optical model analysis is 1083 mb [36]. 
Thus, the CDCC result is consistent with the experimental 
data. Table I shows effects of the spin-orbit interactions and 
the Coulomb breakup on ctr , cteb and cttf- Comparing the 
results with each other, one can find that the Coulomb breakup 
effect is about 50% for (Teb and 2% for ctr, while the spin- 
orbit interaction effect is 17% for cteb and 1% for ctr. Thus, 
the effects are sizable for cteb and appreciable for ctr. Mean- 
while, these effects are quite small for (Ttf, because the abso- 
lute value of S is mainly determined by the imaginary parts of 

U^'' and Up^\ This is also the case with (Tp:STR, ctthStr,, fip 
and (TCP. Table I shows that the spin-orbit interaction effect 
is even smaller than the Coulomb breakup effect. We hence- 
forth neglect the spin-orbit interactions but not the Coulomb 



In this subsection, we discuss the relation between ctr and 
apB for deuteron induced reactions at 200 MeV/nucleon. 

In the framework of CDCC, ctr is the sum of the partial 
reaction cross section (Tb.{J) over the total angular momentum 
J, while (TEB is the sum of the partial breakup cross sections 

crEB(^): 

TR = (•^) = -F2 + (^)' (36) 

,7 J 

TEB = ^'Teb(J) = -^^(2J+1)Peb(J) (37) 
J ^0 J 

with 

Pr(J) = 1-|(0|^(J)|0)|2, (38) 
Peb(J) = ^|(/3|5(J)|0)|2, (39) 

/3 
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where Kq is the initial wave number of d. The partial elas- 
tic and breakup iS-matrix elements are denoted by (0| S'( J) |0) 
and {P\S{J)\0), respectively, where represents the elas- 
tic (breakup) channel. The quantity -Pr( J) shows, for each J, 
the transition probability of the incident flux to all channels 
except the elastic channel, while Peb{J) describes the transi- 
tion probability to all the breakup channels. The probability 
-Feb(^) can be rewritten into 



Peb{J) - {0\S{J)^S{J)\Q) - |(0|5(J)|0)|^ 



(40) 



This indicates that Peb{J) is the fluctuation of the mean 
value |(0|S'(J)|0)| for each J. In general, a rapid change in 
|(0|S'(J)|0)| with respect to J occurs where the fluctuation be- 
comes maximum. Since -Pr(J) is a function of |(0|S'(J)|0)|, 
one can expect that Pr{J) is rapidly changed where Peb{J) 
becomes maximum. We return to this point below. 

The transition probabilities Pr and Peb are plotted in Fig. 3 
as a function of the effective distance R = {J + l/2)/Ko be- 
tween the projectile and the target. For heavier targets, ^^Ni, 
^^Nb and ^"^Pb, Pr behaves as a logistic function and hence 
the R dependence is close to a step function. This indicates 
that the reaction cross section can be approximately described 
by the black-sphere model [37]. Therefore, ctr can be ex- 
pressed by the area of a disk 



(41) 



with effective radius i?R . Meanwhile, the elastic -breakup re- 
action is peripheral, since Peb has a single peak at a finite 
value of R. An effective radius Reb of cteb can be defined 
by the peak of Peb- As expected, Pr changes rapidly at 
R = Reb- This indicates that 



^R — Peb- 



(42) 



For lighter targets such as ^Be and -^^Al, Peb has two peaks; 
the first peak is located at P = and the second at finite value 
of P. However, the second peak is more significant than the 
first peak in ctr because of the weight factor of 2 J + 1 in 
Eq. (36). We thus define Peb by the second peak. 

Figure 4(a) shows Peb as a function of A^^"^, where A is 
the target mass number Since the elastic-breakup reaction is 
peripheral, Peb is expected to depend on A^^^. Actually, A- 
dependence of Peb is well fitted by a straight line (the solid 
line) 



Peb = 0.33 + 1.46Ai/3_ 



(43) 



The fitting is made only for heavier targets of ^^Ni, ^^Nb, and 
^•^^Pb, where Peb{J) has a single peak. Figure 4(b) shows 
A-dependence of <tr , where the solid curve is obtained from 
Eq. (41) with Eqs. (42) and (43), while closed circles stand 
for the results of CDCC. Both results agree with each other 
Thus, the formula (42) is well satisfied. 
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FIG. 3: (Color online) Transition probabilities Pr and Peb as 
a function of i? — (J + 1/2) /Kq for deuteron scattering at 
200 MeV/nucleon. 
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D. A-dependence of integrated cross sections 

A-dependence of integrated cross sections is discussed for 
deuteron induced reactions at 200 MeV/nucleon. Integrated 



TABLE II: Integrated cross sections calculated with CDCC and ERT. 
The cross sections are shown in units of mb. 
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FIG. 4: (Color online) A dependence of (a) 7?eb and (b) ctr. 



cross sections calculated by CDCC and ERT are tabulated in 
Table II. 

First, we consider the total-fusion cross section cttf- The 
cross section is obtained by subtracting the area of the ring 
(Teb (J) from that of the disk ctr ( J). Thus, it can be described 
also by the area of a disk 



(TTF 



ttR^ 



TF 



(44) 



with effective radius Rtf- Similar definition is possible for 

a'TF{p) and <7^p{n): 



c^tf(p) = 7ri?TF(p)^, a-TF{n) = 7ri?TF("-)^ 



(45) 



Figure 5 presents Rtf, Rtf{p) and i?TF("-) as a function 
of yl^/"^. Symbols show effective radii evaluated from the 
CDCC total-fusion cross sections in Table II. They can be 
fitted by straight lines 

Rtf = 0.18 + 1.41^1/3^ (46) 
^tf(p) = -0.60 + 1.36Al/^ (47) 
i?TF(n) = -1.09 + 1.46^1/^ (48) 



This fitting is made only for heavier targets, ^^Ni, ^''Nb, and 
^°®Pb, where Pr{J) has a logistic shape, but the fitting is still 
good for lighter targets of ^Be and ^^Al. 
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FIG. 5: (Color online) yl-dependence of -Rtf, -Rtf(p) and Rt¥(ji)- 
The symbols denote the results of CDCC, while the lines stand for 
the results of the straight-line fitting. 



The neutron stripping cross section is obtained from Rtf 

and Rtf (p) as 

0-„:STR = 7r[i?|p - i?TF(p)^] = 27rAl:STRi?ri:STR (49) 

with 

Dn.STR = Rtf - Rtf (p), (50) 

Jin:STR = [i?TF + i?TF (p)] /2. (51) 

Thus, the neutron stripping reaction occurs on a ring of effec- 
tive radius -R„:str and effective width D„:str - The effective 
width -DniSTR has small A-dependence because of the cancel- 
lation between Rtf and RTpip)- Similar discussion can be 
made for the proton stripping cross section: 



CpiSTR — 27r£)p:STR-Rp:STR (52) 



with 



-Dp:STR = Rtf - RTpin), (53) 

Rp:STK = [i?TF + i?TF(«)]/2. (54) 

The effective radii i?„:STR and Rp.sTB., and the effective 
widths -DniSTR and Dp-sTR are simply obtained from Rtf, 
R'Tf{p), and RTpin): 



i?„:STR = -0.21 + 1.39^1/3, 

i?p:STR = -0.46 + 1.44/11/3, 

^„:STR = 0.78 + 0.05^1/3, 

£'p:STR = 1.27- 0.05A1/3. 



(55) 
(56) 
(57) 
(58) 
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These are shown in Fig. 6 as a function of A^^"^. As expected, 
£'n:STR and Dp-sTR have weak A-dependence. The values 
are about 1 fm corresponding to the diffuseness of the target 
density. Meanwhile, i?,i:STR and i?p:STR has almost the same 
A-dependence as i?R . 
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FIG. 6: (Color online) A dependence of -Rn:STR, J?p;STR, £'„:str, 
and DpiSTR. 



The elastic -breakup reaction is peripheral as the stripping 
reactions. It is thus natural to assume that it occurs on a ring 
with effective radius i?EB and effective width Z^eb^ 



(Tee = 27r£'EB^EB- 



(59) 



The effective width Z3eb may be parameterized by 

Dmi = a + hA^'^ + cZt (60) 

with parameters a, &, and c. Note that the role of Coulomb 
breakup, which is essential for cteb, is described by the last 
term cZt, where Zt is the proton number of target. We use 
the relation between Zt and A for nuclei on stability line: 



Zt 



2+ (ac/aA)A2/3 2 + 0.015^2/3 



(61) 



obtained from the Bethe-Weisacker mass formula [38] ne- 
glecting pairing energy term, where ac = 0.697 MeV and 
OA = 46.58 MeV are coefficients of the Coulomb and asym- 
metry energy terms, respectively. The parameterization (60) 
indeed works well as shown in Fig. 7. The circles denote Z^eb 
evaluated from the CDCC elastic -breakup cross section with 
Eq. (59). A-dependence of the CDCC results is well simulated 
by the solid line 



i:>EB = 0.007 + 0.011^1/3 o.005Zt. 



(62) 



The Deb is found to be smaller than -D„:str and I?p:STR- 
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A 

FIG. 7: (Color online) A dependence of Deb. 



Other integrated cross sections can be obtained by the com- 
bination of o-EB, o-niSTR, crp:STR, and (Ttf: 



(T_„ = 27rZ3EBZ?EB 

CT-p = 27rDEBZ?EB 



27r£)„:STRZ?n:STR, 
27r£)p:STR-Rp:STR,, 



(Tie — 27rD„:STR-Rn:STR + 27r£)p:STRi?p:STR, 

ccF = 7ri?^p — (Tie. 



(63) 
(64) 
(65) 
(66) 



One can thus define effective radii i?cF of the complete fusion 
cross sections by 



O'er 



CF 



(67) 



and can evaluate the value of i?cF from Eq. (66). A- 
dependence of i?TF, Z?n:STR, -RpiSTR, -Reb, and i?cF is 
summarized in Fig. 8. The order of the effective radii is 
i?CF < -Rp:STR, ~ Z?„:STR < i?TF < -Rr = i?EB, inde- 
pendently of A. Among these reactions, the elastic -breakup 
reaction is most peripheral, and it occurs at i?EB — -Deb/2 < 
R < Reb + £'eb/2. The incomplete fusion reactions take 
place at i?„:STR - £',i:Str/2 < R < i?„:STR + Ai:Str/2. 
At i? < RcR, only the complete fusion reaction occurs. 

The cross sections (Jn-.STR, o'p:STR, and tT_p are plot- 
ted as a function of A in Fig. 9. The cross section (t„:str 
((T_„) has similar ^-dependence to (Tp:STR (c-p). The re- 
moval cross sections have stronger ^-dependence than the 
stripping cross sections, since the former include the elastic- 
breakup cross section. 

A-dependence of ctr , cttf, ctif, and ctcf is summarized in 
Fig. 10. For A < 150, fJip is larger than ctcf and becomes the 
largest component of (Tr, while for A > 150, acF becomes 
the largest. 



E. Accuracy of the Glauber model for integrated cross sections 

The accuracy of the Glauber model is investigated for 
deuteron induced reactions at 200 MeV/nucleon. For this pur- 
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FIG. 8: (Color online) A dependence of Reb, Rtf, Rn-.STR, 
Rp-.STR, and RcF. 



pose, we define the relative error 

Sx - [X(CDCC) - X(GL)]/X(CDCC), 



(68) 



where X(CDCC) and X{GL) are integrated cross sections 
calculated with CDCC and the Glauber model, respectively. 
In the Glauber-model calculation, the eikonal and adiabatic 
approximations are made and the Coulomb interaction is set 
to zero. 

Figure 11 shows 5x as a function of A for cteb, cr-n, and 
cr_p in the upper panel and for ctr, cttf, o',i:str,, crp:STR, and 
(TcF in the lower panel. For light targets, say ^Be, the error is 
less than 2% for all integrated cross sections except (tcf- The 
error is 8% for (Jcf, but (Tcf itself is small there. Thus, the 
Glauber model is good at small A, as expected. 

For heavier targets, say ^"*Pb, where the Coulomb breakup 
is essential, the error is 80% for cteb, 20% for cr_„ and a-p, 
and —20% for (Tp:STR- The error is slightly smaller for (T__p 
than for cr_„. However, this is just a result of the cancellation 
in a-p between the positive error for cteb and the negative 
error for crp:STR- The error is less than 6% for (T„:str, cr, 
cttf, and ctcf- Thus, the Glauber model is not good for <teb, 
cr_„, and (T-i,. 



IV. SUMMARY 

The continuum-discretized coupled-channels method 
(CDCC) and the eikonal reaction theory (ERT) are applied to 
d+^^Ni elastic scattering at 200 MeV/nucleon and inclusive 
'^Li{d, n) reaction at 40 MeV. For ctr of the d+^'^Ni scat- 
tering, the CDCC result is consistent with the experimental 
data. The spin-orbit interactions of the proton and neutron 
optical potentials yield a significant effect on the differential 
elastic-scattering cross section, but not on (Tr. For crp:STR of 
the d+^Li scattering, the ERT result is consistent with the 
experimental data. 
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FIG. 9: (Color online) A-dependence of the stripping and removal 
cross sections. (7„:str and crp:STR are plotted in the upper panel, 
while a-n and cr_p are shown in the lower panel. Symbols stand for 
results of CDCC and ERT. Lines denote results of the parameteriza- 
tion (49)-(66). 



^-dependence of several types of integrated cross sections 
is systematically investigated with CDCC and ERT for the 
deuteron induced reactions at 200 MeV/nucleon that corre- 
sponds to typical RIBF and GSI beam energies. The A- 
dependence is clearly explained with simple formulae as fol- 
lows. A black-sphere type reaction such as the reaction, total- 
fusion, and complete-fusion processes occurs on a disk with 
the area of ttRx with effective radius Rx that is well param- 
eterized by ax + 5xj4^/'^. A peripheral reaction such as the 
elastic -breakup, nucleon-stripping, and nucleon-removal pro- 
cesses takes place on a ring 2ttDxRx with effective radius 
and width, Rx and Dx, and A-dependence of Rx is well 
parameterized by ax + bxA^^"^. For neutron- and proton- 
stripping reactions as the incomplete-fusion reaction, the ef- 
fective widths, £'„:STR and DpSTR, are about 1 fm indepen- 
dently of A. The effective radii, i?„:STR and Rp.sTR, are 
smaller than effective radius i?R for the reaction cross sec- 
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FIG. 10: (Color online) A dependence of integrated cross sections. 
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FIG. 1 1 : (Color online) Accuracy of the Glauber model for integrated 
cross sections. 



tion by about 1 fm independently of A, while effective radius 
i?CF for the complete-fusion cross section is smaller than i?R 
by about 2.5 fm independently of A. Thus, ^-dependences 
of i?R, i?n:STR, ^p:STR, and i?cF are simple and similar to 
each other Thus, if ctr, ctcf, CTmSTR, and crp:STR are deter- 
mined experimentally just for two targets, one can estimate 
these cross sections for any target. It is of interest as a future 
work to see whether this property is held for other incident 
energies and projectiles. 

The total-fusion cross section cttf is obtained from mea- 
surable cross sections, ctr and cteb, by ctf = ctr — cteb- 
Similarly, neutron- and proton-stripping cross sections are de- 
termined from measurable neutron- and proton-removal cross 
sections by (t„:str = <^-n - c^eb and (Tp:STR = o-p - cteb- 
It is thus important to determine A-dependence of cteb- How- 
ever, the j4-dependence is known to be complicated [39, 40], 
since it depends on not only A but also the target proton 
number Zt. This problem can be solved by the formula 
(Teb = SttDeb^eb- Effective radius i?EB agrees with 
-Rr = Or + &R A^^"^ with high accuracy, and effective width 
Z^EB is well parameterized by aEB+^EB^^'^'^+CEB^T- Thus, 
A-dependence of cteb is determined, if cteb is measured for 
three targets and ctr is measured for two targets. 

Accuracy of the Glauber model is also tested for the 
deuteron scattering at 200 MeV/nucleon. The accuracy for in- 
tegrated cross sections is summarized as follows. The Glauber 
model is good for light targets, if the interactions between 
projectile and target are clearly determined. For heavy tar- 
gets, however, the model is not good for the elastic -breakup, 
the nucleon-removal, and the proton-stripping cross sections, 
because of the strong Coulomb field, while it is fairly good 
for the other cross sections. It is quite interesting as a future 
work that similar systematic analyses will be made for heav- 
ier projectiles such as Ne and Ca isotopes with larger proton 
numbers. 
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